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ABSTRACT: This contribution deals with identification of fractional-order dy- 
namical systems. We consider systems whose mathematical description is a three- 
member differential equation in which the orders of derivatives can be real numbers. 
We give a discretization method and a numerical solution of differential equations of 
this type. An experimental method of identification is given which is based on eval- 
uation of transfer characteristics. This is a combination of the method of derivatives 
of transfer characteristics and of the method of passive search. The verification was 
performed on systems with known parameters and also on a laboratory object. 

1. INTRODUCTION 

Real objects in general are fractional-order systems, although in some types of sys- 
tems the order is very close to an integer order. Such systems are mainly electronic 
systems composed of quality electronic elements. So far, however, the systems were 
described as integer-order systems, regardless of the negative consequences caused 
by neglecting the real order of the system [4,5]. Disregarding the fractional order of 
the system was caused mainly by the nonexistence of simple mathematical tools for 
the description of such systems. Since major advances have been made in this area 
in the last few years [1,2,3,4,5,6] it is possible to consider also the real order of the 
dynamical systems. Such models are more adequate for the description of dynamical 
systems with distributed parameters than integer-order models with concentrated 
parameters. With regard to this, in the task of identification, it is necessary to 
consider also the fractional order of the dynamical system. In this contribution we 
will concentrate mainly on the identification of parameters (including the order of 
derivatives) for a chosen structure of the model with emphasis on the methods of 
evaluation of transfer characteristics. The verification of the correctness of parame- 
ter identification will be done by using it in systems with known parameters. 

2. DEFINITION OF THE FRACTIONAL-ORDER SYSTEM MODEL 

For the purposes of later development let us consider as the model of a dynamical 
system the following fractional-order three-member differential equation 

a2y^''\t) + my^^Xt) + aoy{t) = u{t) (1) 

where a and P are in general real, with initial conditions y^^\0) = and y{0) = 0. 
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3. NUMERICAL SOLUTION OF THE FRACTIONAL-ORDER 
DIFFERENTIAL EQUATION 



For the numerical calculation of the fractional-order diffrential equation (||) we use 
for an approximation of fractional derivatives in the equation (|l]) the following for- 
mula taken from [6] 

N(t) 

/^\t) ^ it-L)Dty{t) = h-^ b0{t - jh) , (2) 

j=0 



where L is the "memory length", /), is the time step 

N{t) = min 

[z] is the integer part of z, 

J 



bj = i-^yr] (3) 



where is a binomial coefficient. 

The approximation of the equation ([I|) in discrete time steps {i^ = 2,3, ...) 
has the following form [4] 

TO m 

a2h~" bjUm-j + aih~'^ Y ^jVm-j + «01/to = Um (4) 
j=0 j=0 

From the approximation (^) we can derive [4] the following explicit recursive 
formula for the calculation of the values i/m {m = 2, 3, ...) 

^ _ Um - a2/^~" E^l bjym-j - ai/l~^ E^l CjVm-j 

with = 0, yi = 0. An analytical solution of fractional-order differential equations 
is described in detail in [2], and for closed regulation circuits also in [4,5]. 



4. METHOD OF DIFFERENTIATION OF TRANSFER 
CHARACTERISTICS 

Assume the input of the system (|I]) with known values a and (3 is acted upon by 
arbitrary function u{t), and y{t) is the output function of the system, while the 
initial conditions need not be zero in this method. If the input and output functions 
are known functions of time obtained, for example, by measurement, then for given 
discrete times ("^ = 2, 3, ...) we can determine the values of the derivatives of the 
output quantity according to the formula (0). For each time instance tm the equation 
(§) must be satisfied. Therefore after substituting the corresponding values into the 
equation (^) we obtain a system of three linear equations for unknown coefficients 
Oo, tti, By solving this system of equations and obtaining the values of coefficients 
ao,ai,a2 the system, from which we used experimental data in the calculations, is 
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identified. In order to be able to use more measured values Ui, and to make the 
identification more accurate, we consider for the estimate of the vector of unknown 
parameters a as an error criterion the following functional 

T 

E{a) = [ [a2y^"\t) + aiy'^'^\t) + aoy{t) - u{t)]'^dt ^ min (6) 







A necessary condition for the minimum to be achieved is 

dE{a) 



da 



(7) 



Hence, after manipulations we obtain a system of three linear equations for the 
vector of unknown parameters d 



T 

a2 



{y^^\t)fdt + a,jy^''\t)y^^\t)dt + a,j y^^\t)y{t)dt = j y^^\t)u{t)dt 





T T T T 

a2 f y^''\t)y'~^\t)dt + aJ {y^^\t) fdt + ao f y^^\t)y{t)dt = f y^^\t)u{t)dt (8) 



y^''\t)y{t)dt + aJ y^^\t)y{t)dt + aj {y{t)fdt = y{t)u{t)dt 



a2 



where y{t) and u{t) are the measured output and input functions of the system, and 
the derivatives y^"\t) and y^^\t) are the derivatives of the output function of the 
system, obtained by calculation from the formula @. In the discretization of the 
system (H) we substitute the continuous functions y{t),y^"\t), y^^\t), u{t) in discrete 
time instants = 0, 1, 2, 3, M) with their discrete values ym, \ Vm^'^rn and 
the integration we substitute with the summation for the time interval T for which 
we have M + 1 measured values from the system. Then we get 

MM MM 

m=0 m=0 m,=0 m=0 

M M M M 

"2 E y^VI^ + «i E (y'i^f + «o E y^^^ym = E yi^^^m (9) 

m=0 m=0 ■m=0 m=0 

M M M M 

^2 E yin^ym + oi E ym^ym + «o E (ymf = E ymUm 

m=0 m=0 m=0 m=0 

From the system of equations (|^) we can calculate the unknown parameters a, 
and the task of parameter identification of the system is thus solved provided we 
know the values a and 



5. METHOD OF IDENTIFICATION OF FRACTIONAL-ORDER 
DYNAMICAL SYSTEMS 



In general in the sytems of the type ([^) it is necessary to identify the orders of 
derivatives a and (5 which are non-integer and reaL If, under this assumption, we 
used the given method of differentiation of transfer characteristics for the system (|lD, 
we would obtain a system of five nonlinear equations. It is difficult to solve such 
systems of equations, however, therefore we chose a combination of the method 
of differentiation of transfer characteristics for the calculation of the coefficients 
ao,ai,a2 and the method of simple search for the determination of the coefficients 
a and /3. While the first method is analytical, in the second the coefficients are 
calculated iteratively. 

In the identification of the parameters of the model we proceed as follows: 

(a) For the parameters a and /3 we determine the intervals < amm, cimax > and 

< Pmin, Pmax >, in which wc look for the parameters. We divide each of these 
intervals into ria and rijj subintervals 

Ua = (10) 

'^i^Pmax f^min) /t -i\ 

np = (11) 

where e denotes the fineness of division. We choose the starting values a and 
f3 from the chosen intervals. 

(b) We use the method of differentiation of the transfer characteristic to calculate 
for given a and (3 the corresponding values of the parameters Oq, ai, a2. 

(c) We calculate the criterion of the approximation of experimental output values 
of the real system yci and the outputs of the model i/i 

-1 M 

Q = TT-TT.[y^^'y^? (12) 

which serves for the choice of optimum values of the coefficients a and (3. 

(d) If the value of the actually computed approximation criterion is less than the 
value of the criterion from the preceding iteration step, we choose for the 
optimum subintervals those containing last values of the coefficients a and (3. 

(e) After the calculation of the approximation criterion for all subintervals, we 
determine from the optimum subintervals the intervals < (ymin,0(max > and 

< f3min, f3max > and for these subintervals the new subintervals. We repeat 
this procedure until the required accuracy of the calculation of parameters a 
and /3 is achieved, where we have always calculated the values of the parameters 
oo, Qi, 02 by using the method of differentiation of transfer characteristics. 



6. VERIFICATION OF THE METHOD 



The verification of tlie correctness of tlie described metliod of parameter identifica- 
tion will be done by using it in systems with known parameters. 

Consider an integer-order system [4] with coefficients (I2 — I5 CLl — 3, QiQ — 2, 
a = 2, P = 1. The transfer characteristic of this system has aperiodic behavior and 
is shown in Fig.l. For the chosen intervals a G< 2; 2 > and /3 G< 1; 1 > the precise 
values of the coefficients 02 = 1, Oi = 3, Uq = 2 were calculated with the method 
of differentiation of transfer characteristics. For the intervals a G< 1, 5; 2, 55 > and 
P e< 0,7; 1,33 > the following values of the coefficients were calculated with both 
methods as = 1,0001 , ai = 2,99987, oq = 1,99998, a = 1,99993, l3 = 0,99996 , 
which values are very close to the parameters of the identified system. 
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Fig.l Integer-order system 



Fig. 2 Fractional-order system 



Consider a fractional-order system with coefficients [4] 02 = 0,8 , ai = 0,5 , 
uq = 1, a = 2,2 , (3 = 0,9. The transfer characteristic of this system computed 
according to the relation (^), has periodic behavior and is shown in Fig. 2. For 
the chosen intervals a e< 2, 2; 2, 2 > and P G< 0,9; 0,9 > the precise values of 
the coefficients 02 = 0, 8 , ai = 0, 5 , Oq = 1 were computed with the method of 
differentiation of transfer characteristics. For the intervals a G< 1,5; 2, 55 > and 
/3 G< 0, 7; 1, 33 > the following values were computed by using both methods 
a2 = 0, 80005 , ai = 0, 49996 , ao = 0, 99998 , a = 2, 19996 , P = 0, 89989 , which 
values are very close to the parameters of the identified system. 
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Fig. 3 Approximation of system 



Fig. 4 Real object 



Let us consider again the previous fractional-order system with coefficients 02 = 
0,8 , ai = 0,5 , ao = 1, a = 2,2 , /? = 0,9 , whose behavior is given in Fig. 2. 



Let us identify it as a second-order system, which was the classical procedure up 
until now. If we want to identify the system as a second-order system, we have to 
specify the same upper and lower bounds a G< 2; 2 > and /? e< 1; 1 > . With the 
method of differentiation of transfer characteristics the following coefficient values 
were calculated : 02 = 0,76639 , ai = 0,23184 , ao = 1. The system with these 
parameters is actually a classical integer-order model of a fractional-order system; 
their graphical comparison is shown in Fig. 3. Although optically both characteristics 
are rather close to each other, the use of such model in the synthesis of the parameters 
of the regulator is associated with unfavorable consequences for the dynamics of the 
regulation circuit as shown e.g. in [4,5]. 

Let us now identify the parameters of an experimental heating furnace whose 
transfer characteristic is shown in Fig. 4. Assuming it is a second order system 
with a = 2, /3 = 1, we obtain parameters ao = 1,928 , ai = 4892,733 , 02 = 
—73043,36 and the value 1,02.10"^ of the criterion ([12|). For the chosen intervals 
a G< 1,1; 2, 55 > and /3 G< 0,33; 1,3 > , the parameters 02 = —14994,3 , 
tti = 6009,52 , ao = 1,69 , a = 1,31 , (3 = 0,97 were identified for this system, 
and the criterion (|12D had the value 2,7.10"^ (see Fig. 4). As we can see (02) the 
previous two approximations are unsuitable. Assuming a two-member differential 
equation, the parameters take on the values ai = 788,35 , ao = 1,39 , = 0, 73 
and the value of the criterion (0) is 6, 3.10"^. 



7. CONCLUSION 



The above results of parameter identification for fractional- or integer-order dynam- 
ical systems confirm that in known systems the agreement of calculated and actual 
parameters was very good. In the case of real objects identification it is necessary 
to pay attention to the accuracy of the output quantity measurement in the iden- 
tified system, and to the time step of the measurement from the point of view of 
algorithms. 

When identifying the parameters of fractionarl-oder dynamical systems in which 
real orders of derivatives a and (3 are identified, it is suitable, considering the exis- 
tence of local minima of the functional (|^) , to repeat the calculations for different 
intervals < amin,ctmax > and < f3min,Pmax >, in which we look for the parameters, 
and for different fineness of division e in the formulas ([T0|) and (pT]). 
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